March 2001 



Active-sterile neutrino oscillations in the early Universe: 
asymmetry generation at low \6m^\ and the Landau-Zener 

approximation 

P. Di Barii'2 and R. Foot^ 

^ Istituto Nazionale di Fisica Nudeare (INFN), Italy 
^ School of Physics 
Research Centre for High Energy Physics 
The University of Melbourne 
Victoria 3010 Australia 
( dibari,foot@physics.unimelb. edu. au ) 



Abstract 

It is well established that active-sterile neutrino oscillations generate large 
neutrino asymmetries for very small mixing angles (sin^ 29q < 10"^), negative 
values of Sm"^ and provided that \5m?\ ^ 10~'^eV^. By numerically solving 
the quantum kinetic equations, we show that the generation still occurs at 
much lower values of |(^m^|. We also describe the borders of the generation at 
small mixing angles and show how our numerical results can be analytically 
understood within the framework of the Landau-Zener approximation thereby 
extending previous work based on the adiabatic limit. This approximate 
approach leads to a fair description of the MSW dominated regime of the 
neutrino asymmetry evolution and is also able to correctly reproduce its final 
value. We also briefly discuss the impact that neutrino asymmetry generation 
could have on big bang nucleosynthesis, CMBR and relic neutrinos. 
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I. INTRODUCTION 



Previous work has established that z/q, — i^steriie (« = e, /i, r) neutrino oscillations typi- 
cally generate large [0(0. 1)] neutrino asymmetries for negative values of Sm"^ Q and small 
mixing angles (sin^ 2^o ^ 1) Such a large value for the neutrino asymmetry can 

have several interesting consequences and applications including: a suppression of the ster- 
ile neutrino production [^^0, a modification of the standard BBN predictions in different 
ways [3,[f|,^,PHT0[l and, perhaps not emphasized enough, it implies a relic a neutrino back- 
ground in which mainly neutrinos or anti-neutrinos (according to the sign of the generated 
asymmetry) are present. 

We define the asymmetry of the particle species X as the net number of X particles 
at temperature T per number of photons at some fixed initial temperature Tin, such that 

^ Tin ^ ^c\, in- o, comoving volume R^: 

= TJ^ — ' vlJ 

7 

where Nx = nx is the number of X-particles in R^, while is the X-particle density 

The exact borders of the generation in the space of mixing parameters have not yet been 
fully determined. Using the so called "static approximation" in ^ it was found that, to 
have a final neutrino asymmetry larger than 10~^, the following rough lower limit on the 
mixing angle (for a fixed (5m^) holds: 

sin2 2^o^6(5)xlO-iM^j a = e (/z, r) (4) 



^We define (5m^ = — mf , with mi and m2 eigenvalues of the mass eigenstates such that in 
the limit of ~^ they coincide with the a and sterile weak eigenstates respectively. With this 
definition in mind, one can say that negative values of dm^ imply a sterile neutrino lighter than 
the a. 

particular interesting case is the possibility that a neutrino asymmetry generation can solve 
the claimed discrepancy between the value of baryon density from standard BBN with that one 
inferred from recent CMBR anisotropics observations 

^If one defines: 

T.^Tin^^ (2) 

it is simple to chech that the asymmetry Lx is connected to the asymmetry abundance (relative 
to photons) r]x = {nx — ^x)/''^7 by the simple relation: 

r]x = hLx, (3) 

where h = (T^/T)^ < 1 is the dilution factor that takes into account the photon production during 
the electron-positron annihilations around me/2 ~ 0.25MeV. 
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An upper limit on the mixing angle for which the neutrino generation occurs is still un- 
known as numerical solutions are made quite difficult (i.e. CPU time consuming) by the 
presence of a sterile neutrino production prior the onset of neutrino asymmetry generation 
that is strongly coupled to the neutrino asymmetry generation. Moreover, at large angles 
(sin^ 2^0 ~ 10^^), the phenomenon of rapid oscillations of the asymmetry at the onset of 
the generation appears |]12|, making things even more complicated 0. Requiring that the 



final effective number of extra-neutrino species (in energy density) does not violate the BBN 
bound, ANP < AN^^"" < 1, and with the definition /(AA/'^^^) = ln(l - AA"^^^)/ ln(0.4), 
one finds, for a = e (/x, r), the constraint 




sin^ 200 < 2.1 (4.1) x 10"^ —-- f\AN^^''), (5) 



valid for \5m'^\ > 10"^ eV^ (with 5m'^ < 0) and 0.9 > AA"™^^ > 0.3. This constraint 
is however applicable only in the case of two fiavour z/^ — Us (with a = fJ^,T) neutrino 
oscillations. In the general case of multifiavor oscillations the neutrino asymmetry produced 
by the oscillations with the largest |5m^| can suppress the sterile neutrino production from 
the other oscillations Therefore Eq. (^ can only possibly apply for oscillations with 

the largest |5m^|, with the constraints on other two fiavor oscillations being much weaker. 
Also, if a = e, or if the muon or tauon neutrino is also mixed to the electron neutrino then 
an electron neutrino asymmetry can be generated, the constraint AA"^ < AN^^^ has to 
be replaced by a more general constraint AA^^ + AA'/''<= < AN^^^, in which the quantity 
ANI"^ takes into account the BBN effect of the deviation of the electron neutrino distribution 
from the standard case, in which a Fermi-Dirac distribution with zero chemical potential 
is assumed Q. If a positive electron neutrino asymmetry is generated before the freezing of 
the neutron to proton abundance ratio, then AN^"'' is negative |jl|J^J|,|,|ll|] . In this case 
the upper limit on the mixing angle can be strongly relaxed and is related to the limit for 
the neutrino asymmetry generation, that however, as we said, has not been determined. 
Moreover one has conservatively to keep in mind that the current observational situation 
on the values of primordial nuclear abundances cannot strictly exclude a value AA"™'^^ > 1, 
for which no BBN constraints on active-sterile neutrino oscillations would apply at all if 
ANl"^ =0 0. It has also to be mentioned that the region of the mixing parameter space 
in which a positive electron neutrino asymmetry is generated such that AN^"^ is negative 
and with a mixing angle large enough that AA"^ ~ 1, would be very interesting as in this 
case the standard BBN would be modified such that the primordial nuclear abundances 



^Actually the appearance of neutrino asymmetry sign oscillations, in the region where the sterile 
neutrino production is not negligible, may not be a coincidence. It is likely that the strong coupling 
between a copious sterile neutrino production and the neutrino asymmetry evolution is responsible 



for the rapid sign changes [|13 



more detailed discussion on the definition of AA'^^ and AnI"" can be found in pd[ |. 



^More generally no BBN constraint applies if AA'^^'^^ > 1 + max( AA''/"" ) 
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predictions would still agree with the measured values and at the same time a higher value 
for the baryon density would be allowed, in agreement with that one inferred from CMBR 
acoustic peaks |14| , pT| . 
For 

|5m2| > 10- (6) 



the numerical calculations performed with the quantum kinetic equations (QKE) [|15| found 



that the final value of the neutrino asymmetry is in the range 0.23 — 0.35 approximately 
independently of sin 29o for most of the range of values set by Eqs. (^,||) . These calculations 
also showed that the final values can be well reproduced within the adiahatic approximation. 
Let us shortly recall the main features of the neutrino asymmetry generation effect and of 
what the adiabatic approximation consists [|. Neutrino asymmetry evolution in the region 
of mixing parameters set by the equations (|^), (^) and is characterized by the presence 
of a critical temperature that is approximately given by the following expression: 




T, ^Tl(y,) = 15.0 (18.6) ^ - M^V, (7) 



where Uc = Pc/T is the critical (dimensionless) resonant momentum and its value is ~ 2 Q. 
The total effective asymmetry L*^") is defined as 0: 

L(°) =2L,„+L(") (8) 

where L contains the asymmetry of the other neutrino flavours {(3 ^ a) and also a contribu- 
tion coming from electrons and baryons. We assume that L is of the same order of magnitude 
as the baryon asymmetry ~ 5 x 10"^°, a value that we adopted in our calculations The 
evolution of the total asymmetry experiences different stages. At high T ~ Tc, oscillations 



^See also [16| for a more detailed review. 



^This is true until the constraint (^) on the mixing angle is adopted, otherwise yc can grow to 
much higher values [Q,^] as an effect of sterile neutrino production. The asymmetry generation 
should go off when yc becomes so high to be in the tail of the distribution (yc > y"^ ~ 10), but 
this has not yet been studied by numerical calculations and the exact value of y°^ has not been 
determined. 

^When not necessary we will drop the subscript in L^°'\ L^"^ and just indicate them with L,L. 

any case the final value of L is not dependent on a particular choice for L. Moreover note that 
in the realistic case of multiflavor oscillations a mixing I'fjj^a ^ could create an ^ 5 x lO^^'^ 
prior the occurrence of Vs oscillations. Thus, in the case of multiflavor mixing, our assumption 
is valid if either the Va ^ oscillations have the largest |5m^| or the oscillations ^ Vs do not 
have the right parameters (negative 5m^, large enough mixing angle, . . . ) to create a large |-^^^/^|. 
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drive \L\ « \L\. When T approaches a decrease of |L| to arbitrarly small values is pre- 
vented by the presence of chemical potentials in neutrino equilibrium distributions and \L\ 
starts to grow slowly. At T = T^, by definition, L = L and L^^ = 0|^. Until this momentum, 
the resonant momentum ?/res = Pres/T of neutrinos and antineutrinos is practically the same 
and evolving as oc T^^. At T ^ the asymmetry growth becomes exponential and the 
neutrino and antineutrino resonances quickly become separated: for L > the antineutrino 
resonant momentum gets smaller than the critical momentum yc, while the neutrino reso- 
nant momentum increases rapidly moving into the tail of the distribution where the number 
of neutrinos becomes negligibly small (and vice- versa for negative L). When \L^°'^\ reaches a 
threshold value ~ 1.0(0.4) x lO^^ {yc \6m^\/eV'^y/^ for a = e, (/i, r), non linear effects of 
the medium start to act and the exponential growth turns into a power law (with L oc 
initially), while (for L positive) anti- neutrino resonant momentum gets a minimum value 
?/™™ and then starts to grow again. A relation between yj-^s and L can be obtained from the 
resonance condition and is given by the following approximate expression: 



~ a^LiT/MeVy 

with ao = {Ay/2C{?>)/n^) x lO^^^^MeV^ ^ 8. When j/res ~ 10 the asymmetry stops its 
growth and gets frozen to its final value Lf ~ 0.3. It is simple, from the Eq. (|]), to see that 
this happens at a temperature given by: 

T,^0.5Mev(^)'. (10) 

It is quite interesting that this complicated behaviour of the neutrino evolution can be 
regarded as the effect of two different physical regimes. The relevant quantity is the ratio 
r = (£int/^osc)res betweceu the interaction length and the oscillation length for the resonant 
neutrinos. If r ^ 1 collisions are able to inhibit the MSW effect and the regime is collision 
dominated and can be described by a simplified set of kinetic equations obtainable within 
the static approximation procedure p|JT7|. On the other hand if r ^ 1, the MSW effect 



can take place and the regime is MSW dominated. In it has been clarified that the 
transition between the two regimes depends on the value of the asymmetry and not just on 
the temperature. It is in fact possible to show that the following asymptotic expressions for 
r holds: 



sin 26*0 



0.8(2.0)10-2 
sin 29 Q 



0.8(2.0)10-2 \Lt 



m « u) (11) 

(12) 




Within the static approximation a definition of Tc is more rigorous This is a sort of 

numerical definition. 
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First of all note that a collisional regime can exist only if sin 26q ^ 10~^. Thus sin 26o ~ 10~^ 
is a border between two very different regions of mixing parameters. In this paper we will 
only be interested to the case of very small mixing angles (sin 26o -C 10~^) in which collisions 
are important at low values of the asymmetry < Lt), independently on the temperature. 
In this case only when \L\ ^ (10^^/ sin26'o)^''^ the MSW regime can begin. It is then 
crucial that, during the collisional regime, a fast growth of the asymmetry occurs for the 
MSW effect to start to be effective and for the resonance of neutrinos and antineutrinos to 
become separated. If in fact now one assumes (for L > 0) that the expansion is slow enough 
that the resonant momentum of antineutrinos, after the minimum, increases so slowly that 
antineutrinos adiabatically cross the resonance, then all a-antineutrinos are converted into 
antisterile neutrinos. This easily explains why the final values of the a-neutrino asymmetry 
are in a range of values slightly lower than the value Np^(Ti^) / N^(Tin) = 0.375, corresponding 
to an extreme situation in which all anti-neutrinos are converted. This simple picture is the 
adiabatic approximation described in [Q. 

As we said, the numerical solution of the QKE's have shown that for \5m?\ > 10"^ eV^ 
and for the mixing angle within the limits (^) and (^) (without however a systematic inves- 
tigation at the borders) this physical approximated two-regime picture works very well in 
the description of the neutrino asymmetry evolution and of its final value Q. 

What about for \5m?'\ < 10~^eV^ ? Does neutrino asymmetry generation still occur ? 
Does the static- adiabatic approximation picture still work ? 

The authors of ref. |T^, with an analytical procedure based on the mean momentum 



2| > 



approximation |^ were finding that, even though an exponential growth occurs for \Sm 
10~^eV^, soon the non linear term in the effective potential is able stop the growth and 
concluded that no neutrino asymmetry (larger than 10"'') can be generated in any point of 
the parameter space. This result seemed to be confirmed by numerical results [|r^. These 
were obtained also employing the mean momentum approximation for Srn^ = — lO^^eV^ and 
6*0 = 0.1. At some critical temperature an oscillatory behaviour of the neutrino asymmetry, 
with an amplitude exponentially growing, was found, but the amplitude was not growing 
higher than 10~^, although this may simply be because the evolution was stopped at T ~ 
0.3MeV, below the freezing temperature of neutron to proton abundance ratio or simply 
an artifact of not including the momentum degree of freedom. In any case the authors of 
Ref. |T9| concluded that neutrino oscillations cannot generate a large neutrino asymmetry 
in the early Universe, in apparent agreement with UTSf. In |l|,0,^Q it was shown that a large 



neutrino asymmetry can be generated due to the crucial role of the collisions, thus disproving 



^■^Excluding the phenomenon of rapid oscillations at the onset of the exponential growth for 
sin^ 2^0 ~ 10~^ [12| that requires some further physical picture still not completely understood. 

^^This approximation works very well in the description of the sterile neutrino production for 
reasons explained in ||6|. However in the case of the neutrino asymmetry generation, especially in 
the MSW dominated regime, the momentum dependence is such a crucial feature that the mean 
momentum approximation is a very drastic one and in our opinion it should be regarded more as 
a toy model. 
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the general conclusions of [|T8|,|T9[ , although this result is still not necessarily in disagreement 
with the conclusions of [|l|,l£] for \5m?\ ^ 10"*^ and this point has still remained unclear. 

It has also to be said that the main attention was put at higher values of \5m?\ because 
only in this case the asymmetry generation occurs early enough to modify the standard 
BBN predictions for the nuclear abundances. Moreover there are also qualitative reasons 
to suspect that at low \5m'^\ the static approximation is not valid. In it has been 

infact pointed out that for \5m?\ < 10~^eV'^, oscillations betweeen collisions can give a non 
negligible contribution to the rate of the neutrino asymmetry, suggesting a possible deviation 
from the static approximation. This statement is justified observing that, if one requires 
that the change of the mixing angle in matter is negligible between collisions on average 
then the critical temperature has to be higher than approximately 3 MeV and from the Eq. 
(0) this implies just |5m^| ^ 10~'^eV^. On the other hand we already noted that the role of 
collisions in inhibiting the MSW effect, for low values of the asymmetry, is effective at all 
temperatures (for sin 26*0 ^ 10~^). 

Which kind of deviations from the static approximation are then expected at small \6m'^\ 
and, above all, does the asymmetry generation still occur ? Are oscillations in the neutrino 



asymmetry expected to appear as the mean momentum calculations of [|I9| suggest ? 

The purpose of this paper is to answer these questions, mainly understanding the exact 
border of neutrino asymmetry generation at low values of |5m^|. We will however also 
give some results concerning the presence of neutrino asymmetry oscillations. Moreover we 
will study the border of the generation at small mixing angles using the QKE's and show 
that the limit @, derived in the static approximation, has actually to be replaced with 
a less stringent one (section II). We will show that the small mixing angle border can be 
explained within an analysis on the limits of validity of the adiabaticity approximation in 
the description of the MSW dominated regime. This analysis will also involve extending the 
adiabatic approximation in the border region and outside by means of the Landau-Zener 
approximation, that is able to describe how the neutrino asymmetry generation dies at 
small mixing angles (section III). In section IV we will study the cosmological consequences 
of the neutrino asymmetry generation in the light of the new results. If neutrino asymmetry 
generation occurs also at low values of |(5m^|, even though this implies that it is generated 
at temperatures too low to produce remarkable changes in standard BBN, one can still 
consider which kind of cosmological effects (and maybe observations !) can result from 
it. In particular we will correct some wrong claims, which have appeared recently in the 
literature, on the role that the neutrino asymmetry generated in active-sterile neutrino 
oscillations would have on CMBR. We will consider the implications on the possibilities of 
a detection of relic neutrino background by the light of recent studies |^D|]. In section V we 
will draw the conclusions of this work. 



II. NUMERICAL RESULTS 



^^This is equivalent to imposing that the derivative of the effective potential, with negligible 
neutrino asymmetry, is much less that the total interaction rate 
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A. Set of equations and numerical procedure 



We have numerically solved the QKE's for the density matrix distribution in the momen- 
tum space, describing the statistical properties of the z/^, states, including their mixing. 
Since we are interested in extending previous studies to values of |5m^|/eV^ < 10""^ and 
since from the equation ([T(]|), obtained within the adiabatic approximation, one would have 
that the neutrino asymmetry is generated down to temperatures below T ~ 0.05 MeV, we 
have to take into account the occurance of electron-positron annihilations that will make 
Ty 7^ T. The first modification is that the re-scaled momentum pR/Ri^ = pTin/T,^ and thus 
now all neutrino distributions are more conveniently expressed in terms of y = p/T^. 

The QKE's describing neutrino propagation are given by the following set of equations 



for 8 distributions in momentum space |T^, P{y,t), Po{y,t), P{y,t), Po{y,t) (for notational 



simplicity we drop the dependence on momentum and on time): 

dP. 



dt 
dP^ 

dt 
dP^ 

dt 
dPo 

dt 



-XPy-DP, (13) 

\P^-(3P,-DPy (14) 

(3Py + R (15) 

R (16) 



The function R{y,t) takes into account the (thermal) redistribution of neutrinos at each 
quantum state due to the effect of collisions, that in a linear approximation is given by 



R 



Jo 



(17) 



where fp{y,t) = [1 + e^"^^*-*]"^ is the Fermi-Dirac distribution, ^{t) is the (dimensionless) 
chemical potential and /o = /|=o- The variables Pq and Pz can be replaced with the variables 
Zs = {Pq — Pz) 1^2, and Za = {Pq + Pz) /'2- These variables are the diagonal entries of the density 
matrix and are related to the a and sterile neutrino distributions simply by Za,s = fuo,,s/ fo- 
In this case the equations (|T5|) and ( |16|) are replaced by the following ones: 

while the repopulation function R now can be recasted as P = F {z"^ — z^). 

The anti-neutrinos barred distributions P, Pq obey the same equations with the same 
total collision rate F, the same decoherence parameter D, and the same function (3 |^. The 



15 



We neglect a small term, proportional to the asymmetry L, that makes the collision rate F for 
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function A, containing the effective potential of neutrinos, depends explicitly also on L 
and for anti-neutrinos it changes such that 'y{L) = 7(— L). Thus the presence of a total 
asymmetry L couples together the neutrino and antineutrino sets of equations. They are 
coupled also by the presence of the chemical potentials in the equilibrium active neutrino 
distributions. These are infact different if a neutrino asymmetry is present, considering that: 



1 



1 1 



1 + ey-^ 1 + ey~^ 



(20) 



Note that if L = 0, then one would have simply L = 2 L^^ and L = L^^ = would 
be an exact solution of the set of equations. Above the critical temperature this would be 
a stable solution and L would get so small that the statistical fluctuations on very small 
scales would determine the sign of the asymmetry below the critical temperature when the 
solution becomes unstable. This sign would be different on distances as small as the scale 



of the dominant statistical fluctuations that is smaller than the interaction length [|12[. On 
these scales neutrino free-steam and the free-streaming is faster than the neutrino asymmetry 
generation such that immmediately regions with different sign would destroy each other and 
no asymmetry generation could take place. This picture is consistent with the idea that the 
generation of a neutrino asymmetry can occur only if CP is violated at some level, in this 
case due to the presence of a CP asymmetric medium {L ^ 0) that can push the neutrino 
asymmetry toward some direction in any point of the space and strongly enough to dominate 
on statistical fluctuations. This surely happens in a large region of mixing parameters where 
the presence of rapid oscillations can be excluded |T2| . 



Note also that the integral equation (|2(]|) is coupled to the set of differential equations 
for the P, Pq, P, Pq, as the value of the neutrino asymmetry appears in the expression for A. 
However, as already discussed in previous works P,p|,p|,p^ , it is numerically more convenient 



to calculate the neutrino asymmetry not from the equation ( poD but through the following 
integro-differential equation: 



dL,. 



dt 8C(3) 



Jm[Py{y)-Py{y)]My)y'dy. (21) 



The set of equations ( |T^ - |T7| , |2ll) is still not 'closed' as one still needs some equations de- 
scribing the evolution of the chemical potentials ^ and ^ in the equilibrium active neutrino 
distributions. We use the instantaneous chemical decoupling procedure developed in Ref. p| 
with kinetic decoupling temperatures, T^^^ = 3.5 MeV for a = fj,,T and 2.5 MeV for a = e. 
The approximate validity of this procedure has been checked in Ref. 0]. 

The last step is to replace the time with temperature. It is well known that the neutrino 
temperature T^, can be expressed as a function of the photon temperature T, simply imposing 
entropy conservation during electron-positron annihilations: 



neutrinos and f for anti- neutrinos different |^]. Within the static approximation this term proves 
not to give any difference. It has to be investigated whether in the region of rapid oscillations this 
term can play or not any role. 
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gsjmjT) 
9s{0) 



(22) 



with the effective number of (entropy) degrees of freedom of photons and electron-positrons 
given by: 



,2 I ^2 I 1 yl_ 



and 5^5(0) = 11/2 (in figure 1 both and T,y/T are shown for illustrative purposes). One 
can choose either T or Ty as the independent variable for the numerical calculations (we 
used both choices in two different codes). For our discussion it is simpler to choose Ty and 
thus we have to express all quantities as a function of Ty. The quantity (5 can be written as: 

P{y^Ty) = -—sin2e, (24) 

'^y j-u 

Since we treat neutrinos as fully decoupled during the period of electron-positron annihila- 
tions, interactions will be effective only until T ^ T^, and thus we can simply write: 

V{y, Ty) ~ 1.27 (0.92) G^Tly a = e (/i, r) (25) 

The quantity A can be written as: 

\{y, Ty,L) = - —— [cos 26, - b{y, Ty) ± a(y, T,, L)] . (26) 

where the + (— ) sign holds for neutrinos (anti-neutrinos). The term b{y,Ty) is the finite 
temperature contribution to the (dimensionless) effective potential in the early Universe 
pl| . This term is important when \L\ -C Lt, but it becomes negligible after the exponential 
growth, when \L\ ^ Lt, for T ^ Tc. We will mainly consider Tc ^ 0.25 MeV for which 
electron-positron number densities are still approximately equal to the ultra-relativistic limit. 
Thus we can just replace T with Ty and write 0: 

where ^ 18.9 (23.4)MeV for a = e{fi,T). The term a is proportional to the total 
asymmetry L. Electron-positron annihilations do not change electron asymmetry in L and as 
usually is done, we also neglect a small change in L arising from neutron-proton conversions. 



i^Actually we will even consider values of —5m?/eV'^ ^ 10 for which Tc ^ 0.25 MeV. In 
this case one should replace the constant with a function of temperature taking into account 
that electron and positrons number densities decrease during annihilations. We neglected this 
modification and used the standard Notzold-Raffelt expression for the finite temperature term. 
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Thus we can write the term a in the usual way, where now however T has to be replaced 
with Ti, because of the definition (|l]) of particle asymmetry that we are using: 

eV^ 

«^-«o^i^T> (28) 

So far, the electron positron-annihilations have not caused any change in the equations when 
these are written in T^. This because we could approximately assume that a- neutrinos, 
during the electron positron annihilations, are not sensitive to all other particle species, 
either because they do not interact with them any more or because the contribution of the 
other particle species in the effective potentials can be neglected when T < 0.25 MeV and 
\L\ ^ Lt. There is however an indirect effect due to the expansion rate. In replacing the 
time t with the neutrino temperature T,^, one needs to calculate the quantity: 



dt Mpi 



(29) 



where Mpi is the Planck mass and the number of (energy density) degrees of freedom is 
given by []: 

= 2 + -[-) dy -^^^ (30) 

with X = rrie/T and where T has to be regarded as a function of T^, (in figure 1 the function 
Qp is also plotted). Thus electron-positron annihilations cause a change in the expansion 
rate. In section IV we will obtain physical insight into this effect. 

In figures 2,3,4 we show the evolution of the absolute value of the total asymmetry 
with the temperature T. In figure 2 we show solutions for many different values of —Sm^ 
but for a fixed value of the mixing angle sin^ 29o = 10~^ and for a = e. The solid lines 
are the solutions obtained solving the set of QKE equations. The dotted lines are curves 
obtained integrating the QKE's only until r < 3 (see the Eq. |12D. From that moment we 
started to use the adiabatic approximation as described in [0,^. It can be seen how the 
asymmetry generation occurs at values of —5m? much below 10~^ eV^ and goes off only 
around —5m?' ~ 10~^°eV^. In particular for —5m'^ = 10"^'' eV^ the final value of neutrino 
asymmetry is still at the level of 10~^. 

The numerical convergence of the solutions has been obtained integrating around the 
,12] but the interval of the integration highly increases for low \5m'^\ and cov- 



resonance 



ers up to two orders of magnitude in momentum space. The behaviour at low values of 
the asymmetry is well reproduced by the static approximation at any value of |5m^|. In 
particular note that no sign changes are present. 



^^We neglect any contribution to gp from the sterile neutrino production or from the modification 
of a neutrino distribution, due to the generation of the asymmetry. This because we are primarily 
interested in small values of mixing angles satisfying the Eq. (|5|) and to values of \6m?\/eV'^ ^ 100 
for which most of the neutrino asymmetry is generated below T^g^. and thus AiV^ ~ 0. 
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In figure 3 solutions for sin^26'o = 10"''', -Sm'^/eV'^ = 100, 1, 10"^, 10"^ and a = e are 
sliown. Wliile for tlie first two cases tliere are no sign changes, whicli was already known 
1^, in the case —5m? = 10~* two sign changes appear. Note that, changing the initial 



neutrino asymmetry, the solutions converge to the same values at some stage. The final sign 
is still the same of L, like the static approximation predicts. For the numerical convergence 
of this solution a momentum integration over an interval of 5 orders of magnitude has been 
necessary, much larger than for sin^26'o = 10~^. In figure 4 another example is shown 
with the mixing parameters sin^ 2^o = 10^^, Sm^/eV^ = — 10~^ for the a = n case. In 
this case three sign changes are present and the final sign is therefore opposite to that of L. 
We tested this solution with both codes, obtaining the same behaviour. Thus we can say 
that the numerical calculations suggest the appearance of slow oscillations in the neutrino 
asymmetry. These seem to have a different nature than the rapid oscillations observed in 
|ll for sin^ 26^ > 10"^ and 10"^ < \6m^\/eV^ < 500. They are 'slower' and with just a few 
sign changes and do not occur during the exponential growth but before, above the critical 
temperature. Moreover, even though the integration has to be very accurate, a numerical 
convergence can be obtained and thus any character of chaoticity seems to be excluded. 
This is also supported by the fact that starting with two slightly different initial conditions 
the final sign is the same. We can thus exclude the possibility that statistical fluctuations. 



estimated in [12], can influence the solutions. However, it would be desirable to have some 



analytical insight supporting the presence of sign changes found in the numerical results. 
In any case, if confirmed, they would suggest a deviation from the static approximation at 
low \6m^\ and sufficiently high mixing angles (sin^ 

2^0 ~ 10"^) due to the presence of sign 
changes that are not predicted by the static approximation. We did not investigate whether 
at even higher mixing angles, rapid oscillations appear, as in the case of high —5w?. This 
requires a dedicated analysis beyond the goals of this work. 

In figures 5a we give the final value of the asymmetry for a = e, five values of the 
mixing angle (sin^26'o = 10"'', 10"^, 10~^, 4 x 10~^°, 10~^°) and for —5m^ in the range 
10~^^ < —5m?/eV'^ < 10^ 1^. In figure 5b the final value of the asymmetry is shown in 
the case a = fi and for two values of the mixing angle (sin^ 29o = 10~^ and 10~^) and one 
can see that there is not much difference with the case a = e. This because in the MSW 
regime, when the neutrino asymmetry gets its final value, the finite temperature term in the 
effective potential (the 6-term, see Eq. (^) ) is negligible and only the Wolfenstein term is 
relevant (the a-term, see Eq. (^) ). This term is the same for a = e and a = fi. 

For values of the mixing angles sin^26'o = 10~^, 10~^ and 10~^, the final value of the 
neutrino asymmetry can be quite easily understood as being due to MSW transitions as the 
resonance momentum passes through the neutrino distribution. Previous work focussed on 
the adiabatic limit where complete conversion occured. The solid curves describe the results 
obtained solving the QKE. The dotted lines correspond to the results obtained with the 



^^The chosen upper limit for \5m?\ comes from the requirement to have the critical temperature 
Tc ^ 100 MeV considering that, for higher temperatures, the expression that we used for the 
effective potentials are not vahd. However in the case a = e, direct mass measurements from 
Tritium /3 decay hmit —Sm? to be less than ~ 100 eV^. 
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adiabatic approximation, with a slight difference between figure 5a and figure 5b. In figure 
5a {a = e), the results correspond to the solutions in figure 2 and 3 obtained combining 
QKE and the adiabatic approximation as already illustrated. In figure 5b the adiabatic 
approximation is started from T ~ Tc/2 with an initial asymmetry corresponding to a 
resonant momentum y™"^ = 0.3 (see Eq. In the first case one gets a compromise of a 
fast procedure (most of the CPU time is required in the MSW dominated regime) but still 
able to reproduce accurately the correct neutrino asymmetry evolution at any temperature. 
In the second case one has a super-fast procedure able to reproduce accurately the final 
stages (|L| > 10~^) of the evolution, and a fair description of the intermediate stages. At 
high values of — the adiabatic approximation well describes the results obtained using 
the QKE's. Note also that, for \5m?\/eV'^ ^ 10~^, not all anti- neutrinos are converted 
into anti-sterile neutrinos (we refer to positive L, otherwise neutrinos would be converted 
instead of anti-neutrinos) and \L^^\ can be well below the value 0.375 corresponding to a full 
conversion and for ?/™™ -C 1. This is because at high values of |5m^| the MSW dominated 
regime occurs mostly at temperature above ~ 1 MeV and collisions are able to re-distribute 
thermally the asymmetry after it has been produced at y = yj-cs- In particular this means 
that (on statistical average) some anti-neutrinos in quantum states with y > 7/res will move 
to quantum states with y < y^^g. Since y^-es increases with time, this means that these 
anti-neutrinos cannot be converted any more. In a pictorial way, we can say that these 
anti-neutrinos evade the resonance and manage not to be converted. This evasion effect can 
be described correctly only numerically. 

For lower values of —5m? the asymmetry becomes large only for low temperatures ~ 
1 MeV, where collisions have an approximately negligible effect and the evasion effect does 
not occur. The final value of the neutrino asymmetry gets very close to the upper limit 
of 0.375 around —5m? ~ 10~^eV^. For lower values of —5m^ the adiabatic approximation 
predicts that the final value remains approximately constant at 0.375 0. However at low 
enough values of —5m^, the adiabatic approximation clearly fails in reproducing the results 
obtained fully integrating the QKE's. These show that, for a fixed value of sin^ 26'o, there 
is a value of — (5m^/eV^ below which the final value of the asymmetry becomes rapidly 
negligible. The region of mixing parameters where the adiabatic approximation holds and 
the final asymmetry \L^^\ ^0.8 with the value of the adiabatic limit |L|j"|'^^'^^ in 

the range 0.26 — 0.375 (the exact value depending on — 5m^), is given approximately by: 

sin2 2^o f^V~ 1.2x10-9 (31) 



-^^More precisely, as it is shown figure 5a, for sin^ 29q = 10^^ in the adiabatic approximation, the 
final value gets again slowly lower than 0.375 for —5vn?/e\'^ ~ 10^^. This because in that case 
we integrated the QKE's in the collision dominated regime and in this case the correct value for 
the initial y™™ at which the MSW conversions start is reproduced. This value gets as high as 1 at 
—Sw? = 10~^ eV'^ and thus not all-antineutrinos are converted as it happens in figure 5b. On the 
other hand, for larger sin^ 2^0) the value of y'^^ <C 1 and the two procedures, the fast one and the 
super-fast one, give the same answer. 
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It is also simple to infer from figure 5a the following equations describing the iso- 
curves in the space of mixing parameters: 



fin 



sin^2eo ^0.8xlO-^|Lfi:|i 




(32) 



For \L^^\ > 10~^ one finds a larger region compared to the condition (Q) from the static 
approximation. This is because the MSW effect enhances the neutrino asymmetry generation 
compared to a collision dominated regime. The curves (|32D are however valid, if one consider 
values of \L^^\ > 10~^, only for sin^ 29o ^ 10^^. For smaller angles the validity of ( |32D breaks 
down for 10~^ < \Ll'^\ < L*^(sin^ 26^0). For example for sin^ 26^0 = 4 x 10^^° we found that 
L*^ ~ 10^^ corresponding to — 5m^/eV^ = 10~^ and this is clearly visible in figure 5a. For 
lower values of — and a fixed value of sin^ 29q, there is a sharp cutoff in the generation. 

In figure 6 we plotted the iso-|L|^^| curves for — log|L|^^| = 1,2,3,4,5. For sin^ 26*0 ^ 
10^^ the curves (|3^) are found while for smaller mixing angles they do not hold and the 
generation goes off much faster for a decreasing —5rn? and for a fixed value of sin^ 26q. 

One can easily explain this cutoff in the generation of neutrino asymmetry at low values 
of mixing angles (sin^ 26q ^ 10~^) where Eq.'s (0) start not to be valid for \L^\ < L*. This 
can be done in terms of current understanding of neutrino asymmetry generation outlined in 
the introduction. In fact for too small mixing angles the collision dominated regime □ is not 
efficient enough in generating a neutrino asymmetry \L^^^\ ^ (10~^/ sin 26*0)^^^ and the 
MSW dominated regime cannot start at all (see Eq. (|T^) ). This can be seen in figure 7, 
where we plotted the evolution of the neutrino asymmetry with temperature for sin^ 29o = 
4 X 10~^° and a few values of —Sm'^/eV'^. One can see how, in the case — 5m^/eV^ = 10~^, 
the asymmetry gets frozen much before it can get the value (lO'V sin2^o)^/^ - 10~^ and 
the MSW dominated regime never starts. 

We are still lacking in a physical description of the (0) themselves for sin^ 29o ^ 10^^. 
In the next section, we will develop a simple physical picture providing such a description 
improving the adiabatic approximation in the MSW dominated regime. 



We will now develop a useful extension to the adiabatic approximation by taking into 
account non-adiabatic effects which inevitably occur for a choice of mixing parameters out- 
side of the region (|3TD. For definiteness let us refer to positive values of the asymmetry. In 
this case when the total asymmetry has reached values L^"^ ^ Lt, the resonant momentum 



^°We concentrated, for values of the final asymmetry much less than 0.01, on the case a = e. 
However in the case a = /i one does not expect differences, as the Wolfenstein term in the effective 
potential, responsible for the final value of the asymmetry in the MSW regime, is the same. 

^^Note that the rate of neutrino asymmetry generation in the collision dominated regime, as 
deduced within the static approximation, is proportional to sin^ 26*0 



III. THE LANDAU-ZENER APPROXIMATION 
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of anti- neutrinos is given by the expression (^). After it has reached its minimum value y^!^, 
it starts to grow again eventually passing all of the anti-neutrino distribution at y > 
0. In the adiabatic limit all anti-neutrinos at y = y^es are converted. It is then easy to find 
the following rate equation for the neutrino asymmetry 0: 

ad ^ 403) ^^"'^ ~ ^^"'^^ ^ ^^^^ 

In the case of negative asymmetry one has to make the replacement [/p^ — /pj ^ [fus ^ fua\- 
In order to incorporate non adiabatic effects, which can be important if sin^ 26 is small 



enough, we will use the standard Landau-Zener approximation In this case, due to 

the quantum level crossing between matter eigenstates, not all active anti-neutrinos will be 
converted into sterile anti-neutrinos. The number of active anti-neutrinos that undergo level 
crossing and are not converted is given by the well known 'jumping probability' P = Q-'^'y^/'^ 
where 7^ = |2^m^m|rcs is the adiabaticity parameter at the resonance and 9m and im are 
respectively the mixing angle and the oscillation length in matter given by: 

sin 26^ = (34) 
W 2^0 + (cos2^o-&-a)2 



2yT/\Sm^ 



sin2 20o + (cos20o-&-a)^ 



(35) 



Since we are considering a regime of high values of the asymmetry {L ^ Lj), one can 
neglect the finite temperature term b in the effective potential and thus one can find easily 
the following expression for the adiabaticity parameter: 

|(5m2| sin^ 200 , . 

7r = 7^ 7f- I , , ,,| 36) 

2yT.es \da/dt\res 

For computational purposes this can be then transformed with easy algebraic arrangements 
making use of the Eq. (^, defining a = —dlnL^^/dlnTi, > 0, replacing time derivative 
with neutrino temperature derivative (using the expression (^91) ) and finding in the end: 



For a = 4 the adiabaticity parameter is infinite and this is perfectly understandable because 



We recall that at this stage the resonant momentum of neutrinos is by far in the tail of the 
distribution and continues to increase further while the asymmetry grows and therefore neutrinos 
are completely out of the game. 



15 



when this happens, from the Eq. (j^), one can see that ?/res = 1/™™ di/T-es/dt = 0. 
The Landau-Zener approximation can be started from ?/res — Z/™™ ^^is case the rate 

equation for the asymmetry will be given by: 



dt \ dt 



l_e-fT-) (38) 



ad 



This equation can be solved numerically and clearly is much less CPU time consumming 
then the QKE's. The results from the Landau-Zener approximation for the final value of the 
asymmetry are shown in figures 5a, 5b with dashed lines and we find a good agreement with 
the QKE's in both cases a = e and a = fi,T. The Landau-Zener approximation provides 
useful physical insight into the numerical results. First of all we can understand equation 



(pip for the adiabatic region looking at the expression (|37|) for the adiabaticity parameter. 
The jumping probability starts to be non-negligible for 7^ ^ 2. The adiabaticity decreases 
(on average) during neutrino asymmetry (MSW dominated) growth, since initially a = 4 
while when the growth stops at T = Tf one has a = 0. The whole MSW regime will be then 
adiabatic and the final value of neutrino asymmetry will be close to 0.375 if jresiTf) ^ 2. 
Using the expression (0) for the final temperature one can immediately check that this 
condition reproduces the condition (|3lD for the adiabaticity region in the space of mixing 
parameters. We can also try to estimate the equations ( ^2]) for the iso-asymmetry curves. 
For 'y^^ <^ 1 the final neutrino asymmetry is given approximately by: 



Ll: ~ 0.375 I (7.) (39) 

In this equation the quantity (7^.) is the statistical average of the adiabaticity parameter 
done with the Fermi-Dirac distribution with zero chemical potential (assuming that all 
generation of the asymmetry occurs below the chemical decoupling temperature). If one 
neglects the annihilations and assume Ty = T and Qp = 10.75 = const, then this quantity 
can be expressed through the expression (p7| ) where the quantity a must be replaced with 
its averaged value (a) and {L^^) ~ Lf/2. Again, from the expresssion (|) for the resonant 
momentum and using yf^^ ~ 10, obtains an expression for Tf ~ (|5m^|/eV^)^/^ (L|j^)^/^/3, 
and the the iso-asymmetry curves (|3^) are roughly reproduced when (a) ~ 3, that is a 
quite reasonable value considering that a decreases from a = 4 when y = y™™ to a = 
when i/res ~ 10 an the asymmetry gets frozen to its final value. Thus the Landau-Zener 
approximation provides the correct approach to get a physical insight into the generation of 
neutrino asymmetry at low —5m^ and small mixing angles. 



^^Thus in this precise moment, and only in this moment L cx T^^ (exactly). After this moment a 
power-law with a = 4 cannot hold exactly, since if all neutrinos are adiabatically MSW converted, 
the resonant momentum has to start to grow to sustain the asymmetry growth. However while 
L ^ 0.1, L oc (approximately) since the resonance need only move very slowly to produce 
the required asymmetry. One could measure the average value of a fitting the numerical solutions 
if one is really interested in this quantity. See |23| for a detailed analysis on this (in our opinion 
academic) issue. 
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The Landau-Zener approach is also useful to get a physical insight into the effect of the 
electron-positron annihilations. From the espression (^) for the adiabaticity parameter we 
can infact study the following quantity: 



/(T/me) = 



"7r(T = r,)' 




9p 


(1) 


7r 




10.75 





One can plot this function (see figure 1) and discover that it monotonically increases from 1 
to its asymptotical value /(O) ~ 1.1. Thus the answer is that the effect of electron-positron 
annihilatons is to make the MSW conversions less adiabatic and thus to decrease the neutrino 
asymmetry generation. This is an effect clearly negligible in the fully adiabatic region, since 
if 7r(T = Tj,) ^ 1 then also 7^ ^ 1 and the final value of the neutrino asymmetry does 
not change. While, considering the equation (^), one can see that for \Ll2\ < 1 the 
annihilations decrease the final value of 10%, just a correcting effect. 



IV. CAN THE GENERATION OF NEUTRINO ASYMMETRY BE OBSERVED ? 



1. Big Bang Nucleosynthesis 



It remains to briefly discuss the impact of neutrino asymmetry generation on cosmology. 
The most straightforward consequence is a modification of the predictions of standard Big 
Bang Nucleosynthesis |0J^J^. Standard BBN is modified by active - sterile neutrino oscil- 
lations in two ways. One is by changing the energy density, described by AN^, and the 
second is by modifying the standard electron neutrino distributions [e^^'^ + 1]"^, described 
by AN^"^ . This second change occurs both because an electron neutrino asymmetry can 
be generated and also because oscillations below the thermal decoupling temperature at 
~ 1 Mev produce deviations from thermal equilibrium. In the case oi a = fi, r, clearly in the 
idealized case of two neutrino mixing considered in this paper, only the first effect is present 
in first approximation 0. However this effect is not significant unless a significant neutrino 
asymmetry is generated before the chemical decoupling temperature of about 3.5 MeV. This 
in turn requires relatively large values of —5m? ~ 100 eV"^. Even if all the asymmetry is 
produced before the chemical decoupling (for —5m? ^ 100 eV^), one gets at maximum a 
value of AA^^ ~ 0.4 0. In the case a = e the neutrino asymmetry generation gives also a 
significant contribution to AN^"^ and in particular if the generated asymmetry is positive 
then negative values of ANj^"" are possible 0]. However, this effect is only important when 
a significant neutrino asymmetry is produced before the freezing of the neutron to proton 
ratio at T ^ 0.75 MeV. This requires a —5m? ~ 0.01 eV^ |p. This can also be seen from 
figure 5a where the value of the asymmetry at T = 0.75 MeV has been plotted for mixing 



Actually a small effect on electron neutrino distribution, giving a not zero AnI"" , is present also 
for a = /i, r if one describes exactly the chemical decoupling. This because a few electron neutrino- 
antineutrino annihilations will occur to restore /U (or r) number densities depleting electron neutrino 
number densities |24,|]. 
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angles large enough (sin^ 29o > 10~^) that the generation is adiabatic at — 5m^/eV^ ^ 10~^. 
Also in the multi-flavour case where a heavier Uj. (and/or t/^) oscillates with the Ug which 
generates a large L^^ (and/or L^^) asymmetry, and subsequent oscillations of Ur (and/or 
z/^) with Uf, transfers some of this asymmetry to the i>e sector (thereby generating L,^J one 
similarly obtains that significant Lj,^ production requires a -6m^ ~ 0.5 eV^ Thus 
for |5m^| <^ 10^^ eV^, the asymmetry generation effects are unlikely to produce observable 
deviations on the standard BBN predictions. 



2. Cosmic Microwave Background 

It has been pointed out how a neutrino asymmetry generated above the chemical decou- 
pling temperature would imply a AA^^ able to change the matter-radiation equivalence and 
thus modifying the height and position of CMBR acoustic peaks p5|. However in this case 



such an asymmetry cannot be generated from active-sterile neutrino oscillations because, 
as already pointed out discussing BBN, they generate the asymmetry generally below the 
chemical decoupling and thus produce a value of /S.N^ not high enough to produce remark- 
able effects. However a is also produced for mixing angles large enough (see Eq. (|^)) 
that a significant sterile neutrino production occurs. This also modifies the matter-radiation 
equivalence time and gives effects on CMBR anisotropies In the case of negative —5m^ 
the neutrino asymmetry generation can modify the sterile neutrino production, especially in 
multiflavor mixing. Thus this indirect effect of neutrino asymmetry generation has effects on 
CMBR. In other words a calculation of AA^^ including the account of asymmetry generation 
is needed to derive the correct effects on CMBR anisotropies. 



3. Relic Neutrinos 

We have seen that the effects of neutrino asymmetry generation on BBN for low values 
of — are not significant, while on CMBR they can be only indirect. However it is still 
remarkable that neutrino asymmetry generation would lead to a relic neutrino background 
in which mainly anti-neutrinos or neutrinos are present. At the moment nobody has been 
able to experimentally detect relic neutrinos. However it has been noted in pO| that the 



presence of an asymmetry in the relic neutrino background can have observable effects in 
future neutrino detectors. If very high energy a-neutrinos are produced from cosmological 
sources at high redshift, their vacuum probability oscillation can be changed by matter 
effects due to the presence of an asymmetry in the relic background. Thus, if one knows the 
vacuum probability and the fluxes and energy spectrum of produced neutrinos, it is possible 
in principle to observe a change in the oscillation probability induced by matter effects. 
These effects depend crucially on the value of ?7o°^ = Hq L'^°'\ where Hq = 4/11 is the dilution 
factor at the present. Thus the maximum value of ?7o"^ that can be generated from active- 
sterile neutrino oscillations is ~ 0.27. According to the analysis presented in such a 



value can give at maximum effects of about 0.5% in the change of vacuum probability which 
is probably unobservably small. Nevertheless, it shows how the possibilities of detecting 
relic neutrino background can be changed by the presence of an asymmmetry produced by 
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active-sterile neutrino oscillations and thus the issue of relic neutrino background detection 
should be reconsidered in light of this and deserves further investigations. 

V. CONCLUSIONS 

We have shown that the generation of large neutrino asymmetry via active - sterile 
neutrino oscillations occurs for much lower values of \5m^\ than previously supposed. We 
have made a detailed study of the final asymmetry values in the sin^ 29q < 10^^ parameter 
region. The evolution of neutrino asymmetry in this parameter region can be understood 
quite simply by means of the Landau-Zcncr approximation that made possible to derive 
analytical results in good agreement with the numerical ones. For very small —5m? -C 
0.01 eV'^, the oscillations generate the asymmetry too late to have significant effects on BBN. 
The account of neutrino asymmetry generation can modify the sterile neutrino production 
(this is especially true when considered within a multifiavor mixing scheme) and thus has 
to be taken into account when studying the effects on CMBR anisotropics. Properties of 
relic neutrino background are also modified by neutrino asymmetry generation and this 
observation should certainly encourage new investigations on relic neutrinos detection. 
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Figure Captions 

Figure 1. Effect of electron-positron anniliilations on tlie degrees of freedom, on neutrino 
temperature and on the adiabaticity parameter. 

Figure 2. Evolution of L^'^^ = 2 L^^ + L with (photon) temperature for sin^ 29q = 10~^ and 
values of —6m^ as indicated. The solid lines are the QKE's solutions while the dotted lines 
employ an adiabatic description for the growth in the MSW dominated regime. The initial 
value is L(^) = L = 5 x 10-^°. 

Figure 3. The same as in figure 2 but for sin^ 29q — 10~^ and a different choice of values for 
— (5m^. For —Sm^/eV^ = 10~^ we also show an example for a different choice of the initial 
asymmetry. 

Figure 4. An example of evolution with temperature for L^^\ The mixing parameters are 
sin^ 29o = -5myeV^ = 10"^. 

Figure 5. Final values of the neutrino asymmetry in the case a = e (a) and a = (b). 
Solid lines are solutions from the QKE's. Dashed hnes are the solution with Landau-Zener 
approximation. Dotted lines are the adiabatic limit. In the figure 5a we also show with a 
thin solid line the value of neutrino asymmetry at T = 0.75 McV (for mixing angles ^ 10"^) 
that is approximately the freezing temperature of the neutron to proton ratio. 

Figure 6. Adiabatic region and iso-|L|^°| curves. 

Figure 7. The same as in figure 2 but for sin^ 29o = 4 x 10~^° and a different choice of 
values for —5m?'. For —Sm^ = 10~^eV^ the MSW dominated regime never starts. 
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figure 4 
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